RKKY interaction in carbon nanotubes and graphene nanoribbons 
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We study Rudermann-Kittel-Kasuya-Yosida (RKKY) interaction in carbon nanotubes (CNTs) 
and graphene nanoribbons in the presence of spin orbit interactions and magnetic fields. For this we 
evaluate the static spin susceptibility tensor in real space in various regimes at zero temperature. In 
metallic CNTs the RKKY interaction depends strongly on the sublattice and, at the Dirac point, is 
purely ferromagnetic (antiferromagnetic) for the localized spins on the same (different) sublattice, 
whereas in semiconducting CNTs the spin susceptibility depends only weakly on the sublattice and is 
dominantly ferromagnetic. The spin orbit interactions break the SU(2) spin symmetry of the system, 
leading to an anisotropic RKKY interaction of Ising and Moryia-Dzyaloshinsky form, besides the 
usual isotropic Heisenberg interaction. All these RKKY terms can be made of comparable magnitude 
by tuning the Fermi level close to the gap induced by the spin orbit interaction. We further calculate 
the spin susceptibility also at finite frequencies and thereby obtain the spin noise in real space via 
the fluctuation-dissipation theorem. 

PACS numbers: 73.63. Fg, 71.70.Gm, 75.30.Et, 75.30.Hx 
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I. INTRODUCTION 

The Rudermann-Kittel-Kasuya-Yosida (RKKY) inter- 
action is an indirect exchange interaction between two 
localized spins induced by itinerant electrons in a host 
material^— This effective spin interaction, being deter- 
mined by the static spin susceptibility, is not only a fun- 
damental characteristics of the host system but also finds 
interesting and useful applications. One of them is the 
long-range coupling of spins between distant quantum 
dots^£ which is needed in scalable quantum computing 
architectures such as the surface code^ built from spin 
qubitSfi In addition, the RKKY interaction, enhanced 
by electron-electron interactions, can initiate a nuclear 
spin ordering that leads to striking effects such as he- 
lical nuclear magnetism at low temperatures Such a 
rotating magnetic field, equivalent to the presence of a 
uniform magnetic field and Rashba spin orbit interaction 
(SOI) in one-dimensional systems,— is interesting for Ma- 
jor ana fermion physics in its own right The RKKY 
interaction, proposed long ago for normal metals of Fermi 
liquid type,— ~— ^ was later extended in various ways, in 
particular to low-dimensional systems with Rashba SOI 
in the clear>i£ and the disordered^ limit, and to systems 
with electron-electron interactions in one^ii and two 
dimensions within and withou t 19 ' 20 Rashba SOI. Also, 
a general theorem of Mermin- Wagner type was recently 
proven for isotropic RKKY systems that excludes mag- 
netic ordering in one and two dimensions at any finite 
temperature but allows it in the presence of SO L 21 ' 22 

Moreover, due to recent progress in magnetic nanoscale 
imaging^ one can expect that the direct measurement 
of the static spin susceptibility has become within exper- 
imental reachj 2 ^ For this latter purpose, graphene offers 
the unique advantage over other materials such as GaAs 
heterostructures in that its surface can be accessed di- 
rectly on a atomistic scale by the sensing device. All this 
makes the spin susceptibility and the RKKY interaction 



important quantities to study. 

Recently, the RKKY interaction in graphene has at- 
tracted considerable attention! 2 ^— Graphene is known 
for its Dirac-like spectrum with a linear dispersion at low 
energies. This linearity, however, can give rise to diver- 
gences in the expression for the spin susceptibility in mo- 
mentum spaced and complicates the analysis compared 
to systems with quadratic dispersion. However, Kogan 
recently showed that these divergences can be avoided by 
working in the Matsubara formalism.— 

In the present work we consider the close relatives 
of graphene^i namely carbon nanotubes (CNTs) and 
graphene nanoribbons (GNRs), with a focus on spin or- 
bit interaction and non-uniform magnetic fields. Metallic 
CNTs also have a linear spectrum, so for them the imagi- 
nary time approach developed in Ref. [3^ is also most con- 
venient and will be used here. Analogously to graphene, 
we find that the static spin susceptibility changes sign, 
depending on whether the localized spins belong to the 
same sublattice or to different sublattices^ No such de- 
pendence occurs for CNTs in the semiconducting regime, 
characterized by a gap and parabolic spectrum at low fill- 
ings. 

The spin orbit interaction in CNTs is strongly en- 
hanced by curvature effects in comparison to flat 
graphenef^r— while in GNRs strong SOI-like effects can 
be generated by magnetic fields that oscillate or rotate in 
real spacej 12 ' 13 Such non-uniform fields can be produced 
for instance by periodically arranged nanomagnctSj2& 
Spin orbit effects break the SU(2) spin-symmetry of 
the itinerant carriers and thus lead, besides the effec- 
tive Heisenberg interaction, to anisotropic RKKY terms 
of Moryia-Dzyaloshinsky and of Ising form. Quite re- 
markably, when the Fermi level is tuned close to the 
gap opened by the SOI, we find that the isotropic and 
anisotropic terms become of comparable size. This has 
far reaching consequences for ordering in Kondo lattices 
with RKKY interaction, since this opens up the possibil- 
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ity to have magnetic phase transitions in low-dimensional 
systems at finite temperature that are tunable by electric 
gates. 

We mention that similar anisotropics have been found 
before for semiconductors with parabolic spectrum and 
with Rashba SOI in the clean-^ and in the disordered^ 
limit. However, the spin orbit interactions in CNTs and 
in GNRs are of different symmetry and thus both of these 
problems require a separate study, apart from the fact 
that the spectrum is linear. 

For all itinerant regimes we consider, the RKKY inter- 
action is found to decay as 1/R, where R is the distance 
between the localized spins, thus following the standard 
behavior for RKKY interaction in non-interacting one- 
dimensional systems.— [In interacting systems, described 
by Luttinger liquids, the decay becomes slower] 9 ' 17 ' 37 ] In 
contrast, the overall sign as well as the spatial oscillation 
periods of the RKKY interaction are non-generic and de- 
pend strongly on the system and the regimes considered. 

Finally, we will also consider the dynamical spin sus- 
ceptibility at finite frequency. Via the fluctuation- 
dissipation theorem we obtain from this the spin- 
dependent dynamical structure factor in position space, 
which describes the equilibrium correlations of two local- 
ized spins separated by a distance R. 

The paper is organized as follows. Sec. [TT] contains 
different approaches to the RKKY interaction includ- 
ing imaginary time formalism for metallic CNTs with 
a linear spectrum and the retarded Green functions in 
the real space formalism for semiconducting CNTs. In 
Sec. IIIII the low energy spectrum of CNTs is shortly 
discussed. Afterwards the spin susceptibility is calcu- 
lated both in the absence of the SOI (Sec. IIV[) and in 
the presence of the SOI (Sec. [Vj. In addition, in Sec. 
I VII we present results for the case of a magnetic field 
along the nanotube axis. Such a field breaks both orbital 
and spin degeneracy, leading to non-trivial dependence 
of the spin susceptibility on the chemical potential. The 
fluctuation-dissipation theorem connects the spin suscep- 
tibility and the spin fluctuations, allowing us to explore 
the frequency dependence of the spin noise at zero tem- 
perature in Sec. IVIII The RKKY interaction in armchair 
graphene nanoribbons is briefly considered in Sec. IVIIII 
Finally, we conclude with Sec. IIXI in which we shortly 
summarize our main results. 



II. FORMALISM FOR RKKY 

The RKKY interaction 1 " 3 was studied for a long time 
and several approaches were developed. In this section 
we briefly review those used in this work. 

The RKKY interaction is an effective exchange inter- 
action between two magnetic spins, Ij and lj, localized at 
lattice sites Ri and Rj, respectively, that are embedded 
in a system of itinerant electrons with spin- 1/2. These 
electrons have a local spin-interaction with the localized 



spins, described by the Hamiltonian 

H int = jJ2 s r I h (1) 

l=i,j 

where Hsi/2 is the electron spin operator at site R/, and 
J is the coupling strength. Using second order perturba- 
tion expansion in J^ i 14 ' 15 ' 30 the RKKY Hamiltonian^ — 
becomes 

Hrkky = J 2 X%Ifll (2) 

where x^p = Xaffi^Uj,^ = 0) is the static (zero- 
frequency) spin susceptibility tensor, and where summa- 
tion is implied over repeated spin indices a, j3 = x,y,z 
(but not over Here, we assumed that the system 

is translationally invariant so that the susceptibility de- 
pends only on the relative distance Ry = R^ — Hj . The 
RKKY interaction can be expressed in several equiva- 
lent ways. For example, in terms of the retarded Green 
function G(Ry ■: e + i0 + ) the RKKY Hamiltonian is given 
by 

Hrkky = -— Im f'&Tr [(I 4 • s)(I J • s) 

xG(R ij ;e + iO + )G{-R ij ;e + iO + )] , (3) 

where the integration over energy e is limited by the 
Fermi energy tp (see Ref. [Tol ). Here s is a vector of 
the Pauli matrices acting on the spin of the itinerant 
electrons, and the trace Tr runs over the electron spin. 
The retarded Green function G(R;_,-; e + i0 + ), which are 
spin-dependent here and represented as 2 x 2-matrices in 
spin space, are taken in real (Ry) and energy space (e). 
In the presence of spin orbit interaction, we will use Eq. 
^ as a starting point. 

In the absence of spin orbit interaction, the spin 
is a good quantum number, so the effective Hamilto- 
nian Hrkky can be significantly simplified, x*£p = 

6 a pXo i an d the RKKY interaction is of Heisenberg type 
(isotropic in spin space). Expressing the Green func- 
tions in terms of the eigenfunctions ip n (Ri) of the elec- 
tron Hamiltonian, we obtain 

X« = 2 ^ (RiJVm (R< W« (R, Wm (R* ) 

x n F (£ n ) - n F (£ m ) ^ 

where the sum runs over all eigenstates of the spinlcss 
system, and the factor 2 accounts for the spin degeneracy. 
The energy is calculated from the Fermi level, = e n — 
€f, and the Fermi distribution function at T = is given 
by n F (0 = 0(-0. 

The sum in Eq. ([4J is divergent in case of a linear 
spectrum i 30 ' 38 To avoid these divergences we follow Ref. 
|30| and work in the imaginary time formalism, again ne- 
glecting the spin structure of the Green functions. The 
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static real space spin susceptibility at zero temperature 
is given by 



in 

Xo 



2 f°° 

-J dr Go(Rij)T")Go(-R»j,-r) 



(5) 



where the factor 2 again accounts for the spin degeneracy. 
The Matsubara Green functions for t > are found as 

G (R !; „±t) =^J2^n( R ^CR j )e^ T/h e(±^ n ). (6) 



All three approaches to the RKKY interaction described 
above [see Eqs. (J3j> , (jU), and fl5j] arc equivalent to each 
other. Which one is used for a particular case depends 
on calculational convenience. 



where v — 1 (y = —1) corresponds to electrons (holes), 
and n = (k, 7, v) labels the eigenstates. The correspond- 
ing wavefunctions with sublattice spinor arc given by 



i( 7 K+k) 



k G + ik 



R; _J_ f 1 

V2 W e * 



k 2 



(9) 
(10) 



where k = (k G , k). From now on we redefine Dirac points 
by shifting the circumferential value of K (K') by kc, 
so that K = — K' = 47r(tcos# + zsin#)/3a + tk G and 
k= (0,/fe). 



B. Spin orbit interaction 



III. CARBON NANOTUBES 

In this section, we discuss the effective Hamiltonian 
for a carbon nanotube. A carbon nanotube is a rolled-up 
sheet of graphene, a honeycomb lattice composed of two 
types of non-equivalent atoms A and B. The (N±,N2)- 
CNTs can be alternatively characterized by the chiral an- 
gle 9 and the diameter d— The low-energy physics takes 
place in two valleys K and K'. These two Dirac points 
are determined by K = — K' = 47r(tcos# + zsin#)/3a, 
where a is the lattice constant. The unit vector z points 
along the CNT axis, and t is the unit vector in the trans- 
verse direction. 



A. Effective Hamiltonian 

In the absence of spin orbit interaction CNTs arc de- 
scribed by the effective Hamiltonian Hq, 

H = hv F (k G ai + kja 2 ). (7) 

The Pauli matrices ct 2 ; act in the space defined by the 
sub-lattices A and B. The Fermi velocity in graphene 
vp is equal to 10 6 m/s. Here, 7 = 1 (7 = —1) labels 
the K (K') Dirac points, and k is the momentum along 
the z-axis calculated from the corresponding Dirac point . 
The momentum in the circumferential direction k G is 
quantized, k G — 2(m — jS/3)/d, with d the CNT di- 
ameter, leading to two kinds of nanotubes: metallic and 
semiconducting. Here, m £ Z is the subband index and 
6 = (iVi - N 2 ) mod 3 for a (N l7 A^ 2 )-CNT (see Ref. SH). 
The spectrum of metallic CNTs (with k G — 0) is a Dirac 
cone, i.e. linear and gaplcss. In contrast to that, the 
spectrum of semiconducting CNTs (with k G 7^ 0) has a 
gap given by 2hvp\kQ\. In the following we consider only 
the lowest subband with energies 



e„ = vhv F Jk 2 + k%, (8) 



Spin orbit interaction in nanotubes arises mostly from 
curvature effects, which substantially increase its value 
in comparison with flat graphene.— ~ 35 i 39 ~ — The effective 
Hamiltonian, which includes the spin orbit interaction 
terms, is given by 

H so = H + aais z + ~fj3s z , (11) 

where Sj are the Pauli matrices acting on the spin. The 
SOI is described by two parameters, a and ft, which 
depend on the diameter d. The values of these pa- 
rameters can be found in the framework of the tight- 
binding model, a = —0.16 meV / d [nm] and ft — 
-0.62 mcVcos(36>)/d [nm] (see Refs. H3 and SI) . The 
valley index 7 and the spin projection on the nanotube 
axis s are good quantum numbers due to the rotation 
invariance of the CNT. The conduction band spectrum 
{y = 1) is given by 

e„ = ±7/3s + ^J(hv F k) 2 + (~/hv F k G + as) 2 , (12) 
and the corresponding wavefunctions are given by 

^(R l ) = e l(7K+k) - R '-^( et i„)k>, (13) 

e i+.„ = 7fcc + sa + ijk 

^k 2 + ( 7 fc G + so) 2 ' 

where the index n = (k, 7, s) labels the eigenstates. Here 
|s), the eigenstate of the Pauli matrix s z , corresponds to 
the spin state with spin up (s = 1) or down (s = — 1). 
We note that the SOI lifts the spin degeneracy and opens 
gaps at zero momentum, k = 0. In case of semiconduct- 
ing nanotubes (k G ^> k), we use the parabolic approxi- 
mation of the spectrum, 

e„ = hv F k G + -fs(ft + a) + hv F k 2 /2k G . (15) 

Further we denote the sum of the SOI parameters a 
and ft as ft + = ft + a. Such kind of a spectrum is 
similar to the spectrum of a CNT in the presence of a 
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pseudo-magnetic field that has opposite signs at oppo- 
site valleys. We note that there is a principal difference 
between a semiconducting CNT and a semiconducting 
nanowire with Rashba SOI. In the latter, the Rashba 
SOI can be gauged away by a spin-dependent unitary 
transformation . 10 i 15 In contrast, the spectrum of CNTs 
consists of parabolas shifted along the energy axis and 
not along the momentum axis as in the case of semicon- 
ducting nanowires, so the SOI cannot be gauged away. 
As shown below, this leads to a less transparent depen- 
dence of the spin susceptibility on the SOI compared to 
nanowires 



IV. RKKY IN THE ABSENCE OF SOI 
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FIG. 1. The Dirac spectrum of a metallic nanotube. Each 
branch is characterized by the isospin value a, which is an 
eigenvalue of the Pauli matrix 02 ■ We note that the Kramers 
partners at K and K' , respectively, are characterized by the 
same value of the isospin, and two partner states at the same 
cone are characterized by opposite isospins. 



In this section we calculate the spin susceptibility ne- 
glecting spin orbit interaction, so all states are two-fold 
degenerate in spin. We can thus consider a spinless sys- 
tem and account for the spin degeneracy just by intro- 
ducing a factor of 2 in the expressions for the spin sus- 
ceptibility, see Eqs. (01 and ([5]). 



A. Metallic nanotubes 

The spectrum of a metallic nanotube is linear, see Eq. 
(|5|), with the momentum in the circumferential direction 
kc equal to zero, kc = 0. As was mentioned above, in 
this case the integrals over the momentum in Eq. (j4|) are 
divergent^ so it is more convenient to work in the imag- 
inary time formalism [see Eq. ([5])], where all integrals re- 
main well-behaved. To simplify notations, we denote the 
distance between the localized spins as R. — R^ — R^ - 
and its projection on the CNT axis as z. Using the 
wavefunctions given by Eq. Q, we find the Green func- 
tions from Eq. where we replaced sums by integrals, 
J2k ( a /2 7r ) / dk. The Green functions on the same 
sublattices are given by 



straightforward integration 



Gtf A (R, t) = (R, r) = — cos(K • R) 
vpr cos(kpz) — zsa\(kpz) 



(16) 



(vft) 2 + z 2 

The Green function on different sublattices is given by 

G^ s (R,t) =i -sin(K-R) 

Vft sin(kpz) + z cos(kpz) , 
X 'l V2 i 2 • ' > 

Here, the Fermi wavevector kp is determined by the 
Fermi energy tp as kp = tp/hvp. For the corresponding 
spin susceptibilities [see Eq. ((SJ] we then obtain after 



x AA (R) 



4:irhvp\z 



Airhvp\z 



[1 +cos(2K • R)]cos(2/c F z), (18) 



[l-cos(2K-R)]cos(2/c F z). (19) 



If the chemical potential is tuned strictly to the Dirac 
point, kp = 0, the spin susceptibility is purely ferromag- 
netic for the atoms belonging to the same sublattices, 
Xo-iXo B < 0j whereas it is purely antiferromagnetic for 
the atoms belonging to different sublattices, \o B > 0^ 
For the chemical potential tuned away from the Dirac 
point we observe in addition to the sign difference oscil- 
lations of the spin susceptibility in real space with period 
of half the Fermi wavelength ir/kp. This oscillation, to- 
gether with the l/z decay, is typical for RKKY interac- 
tion in one-dimensional systems 

An immediate consequence of the opposite signs of the 
susceptibilities in Eqs. (|19[) is that any ordering of spins 
localized at the honeycomb lattice sites will be antifer- 
romagnetic. Such order produces a staggered magnetic 
field that can act back on the electron system and give 
rise to scattering of electrons between branches of oppo- 
site isospin a at the same Dirac cone (see Fig. [TJ. It 
has been shown elsewhere that such backaction effects 
can lead to a spin-dependent Peierls gap in the electron 
system.— _ — 



B. Semiconducting nanotubes 

1. Zero chemical potential 

Now we consider semiconducting CNTs that are char- 
acterized by a non-zero circumferential wavevector kc 
and a corresponding gap in the spectrum. We begin with 
the case of the Fermi level tp lying in the middle of the 
gap, tp — 0. As a result, there are no states at the Fermi 
level. This leads to a strong suppression of the RKKY in- 
teraction. For example, the Green function on the same 
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sublattice, found from Eqs. (|6]) and ([9]), is given by 



G aa (R,t) 



2tt 



cos(K • R) e «tt*-«fMka 



/27rfc 



£ - z 2 fe G /2-i;j,|r| 
1 



(20) 



where we used the simplified parabolic spectrum, see Eq. 
(If 51) . The spin susceptibility is then obtained from Eq. 



x AA (R) 



2a 2 fc G 
ttHvf 



cos 2 (K-R)A" (2%/2fc G |z|), (21) 



where A'o(a;) is the modified Bessel function of sec- 
ond kind, whic h deca ys exponentially at large distances, 
Kq(x) rs e~ x y / 7r/2x for a; ^> 1. This exponential decay 
(on the scale of the CNT diameter d) of the spin suscep- 
tibility in the case when the Fermi level lies in the gap 
is not surprising. There are just no delocalized electron 
states that can assist the effective coupling between two 
separated localized spins. From now on we neglect any 
contributions coming from higher or lower bands. 



2. Non-zero chemical potential 

In this subsection we assume that the Fermi level is 
tuned in such a way that it crosses, for example, the 
conduction band. Thus, the itinerant states assist the 
RKKY interaction between localized spins. The spec- 
trum of a semiconducting CNT is parabolic [see Eq. 
(|15l) ]. so it is more convenient to work with the spin 
susceptibility given by Eq. ((4]). In momentum space 
the spin susceptibility on the same sublattice is given by 



• M (9)=E 1 , f X M (7,l';9) with 



,AA 



(7>7';g) = 9 1] 



M fc + «)-6r(*) 



(22) 



Performing integration over momentum k, we arrive at 
an expression similar to the Lindhard function, 



x aa (i,i';q) 



akr 



2-nhvpq 



In 



2k f 



2k f 



(23) 



where the Fermi momentum kp is defined as kp = 
^/2fcc(ei? — hvpko)jhvp. Next we go to real space by 
taking the Fourier transform of X AA {lt 7'i ?)• This can 
be readily done by closing the integration contour in the 
upper (lower) complex plane for z > (z < 0) and de- 
forming it around the two branch cuts of the logarithm 
in Eq. (f2"51) that run from ±2kp to ±oo. This yields, 



x AA (K) = 



a 2 kg 



si(2k F \z\) [1 + cos(2K • R)] . (24) 



nhvp 

Here, the sine integral is defined as 

/ smr 7r 
*{x)= dt—--, 



(25) 



and at large distances, x^> 1, its asymptotics is given by 
si(x) ~ — cos(x)/x. 

The evaluation of the spin susceptibility for different 
sublattices is more involved, 



x AB (7,7';?)^£' 

k 



-U* h q O[^(k + q)]-0Mk)] 

£r(fc)-M* + 9) ' 
(26) 

where the phase difference given by A(f)k, q = 4>k — 4>k+q 
depends on the momenta k and q, see Eq. (|TU)) . Tak- 
ing into account that kc is the largest momentum 
characterizing the system, we expand the phase factor 
as e l ^ k ~^ k+q ^ ~ 1 — iq/kc- The main contribution 
to the spin susceptibility comes from the momentum- 
independent part and is given by 



x AB (R) 



a 2 kg 

Tthvp 



si(2fc F |z|)[l-cos(2K-R)]. (27) 



In the next step we evaluate the correction Ax AB (R) 
to the spin susceptibility x AB (R) arising from the 
momentum-dependent part in the phase factor e lA ^ k - q . 
In momentum space it is given by 



2irhvp 



In 



2ki 



q-2ki 



(28) 



By taking the Fourier transform of Eq. 
at the following expression, 



we arrive 



A X AB (R) 



7,7' 



.( 7 -7')K-R > 2 77' 

2irhvp 
<1 - 



dq e lqz In 



2k t 



<1 



2k T 



(29) 



The integral in Eq. (|29|) can be evaluated easily by rec- 
ognizing it as the derivative of si(a;) [see Eq. ([23ll ]. 



/ = 



dx 



In 



1 



. sin a 



(30) 



where a is real. As a result, the correction to the spin 
susceptibility on different sublattices x AB is given by 



A X AB (H) 



a 2 sin(2kpz) 
irhvp\z\ 



[1 - cos(2K • R)] 



(31) 



We note that Ax AB (R.) is small in comparison with 
X^ S (R) by a factor kp/kg *C 1, and, thus, this cor- 
rection plays a role only around the points where the os- 
cillating function ^^(R) vanishes. We emphasize that 
the spin susceptibility for semiconducting CNTs does not 
possess any significant dependence on the sublattices in 
contrast to metallic CNTs. 



V. RKKY IN THE PRESENCE OF SOI 

In the presence of spin orbit interaction, the spin space 
is no longer invariant under rotations, and as a conse- 
quence the spin susceptibility is described by the tensor 
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FIG. 2. The spin susceptibility Xa/3, given in Eqs. (|44[) - (|46p . plotted as a function of distance between localized spins z, for a 
semiconducting (11,1)-CNT in the presence of spin orbit interaction and at zero B-field. Here, \o = a 2 kG /2nhvF ■ For clarity 
we suppress the fast oscillating factors and plot only the slowly varying envelopes at small (a,c) and large scales (b,d). The 
chemical potential [i — 472 meV (<5/x = [i — hvFkc + /3+ = 1 meV) is tuned inside the gap opened by SOI with corresponding 
value 2/3+ = 1.7 meV for a (11,1)-CNT.— The diagonal components \ zz [(a,b) blue full curve] and Xn ( c id) oscillate with 
period 2ir/k+ t - ~ 2n/kF- In contrast to that, the off-diagonal component \xy [(a^b) red dashed curve], oscillates with period 
2n/k+,+ ~ it/ltF- All components decay as 1/z. Note that the diagonal and off-diagonal components are of comparable size. 



Xa/3 [ see ^q- ©] with non-vanishing off-diagonal com- 
ponents. In this case it is more convenient to work in the 
framework of retarded Green functions^ in which the 
RKKY Hamiltonian Hrkky is given by Eq. Be- 
low we neglect the weak dependence of the susceptibility 
on the sublattice discussed above and focus on the SOI 
effects. The Green functions in the energy-momentum 
space can be expressed as 

G(k, e + i0+) = 5}G (7, k, e) + Gi( 7 , k, e)s z ], (32) 

7 

where the diagonal and off-diagonal (in spin space) Green 
functions are given by 



r/— ±1 



frv F ^ fc2 - fc2 + i0 + ' 

7] — ±1 



(33) 
(34) 



Here, to simplify notations, we introduced the wavevector 
k n i, defined as (with ?/ = 777) 



Ay = y/2k G (e - hv F k G - Ti'/3+)/hv F , (35) 

which can take both real and imaginary values. In a next 
step we transform the Green functions from momentum 
to real space, 

GV(R,e + iO+) = |X; / ^oafT.Mle^, 
7 ?r 

(36) 



leading to 
G (R,e) 

Gi(Re) 



. ak G 
-1 — - — cos 

2nvF 



(K • R) 



ak G ■ l\r -r>\ ^ 

sm(K ■ R) > 

2hv F K ' ^ k„ + i0+ 



(37) 



(38) 



=±1 ""J 



Substituting G(R, e) = G (R, e) + Gi(R, e)s z into Eq. 
([3]), we find for the RKKY Hamiltonian, 



J 



Hrkky — Im 

7T 



4/fJ/ f'de G?(R,e) 

J — OO 

4i(I(Xl,-)x / G 1 (R,e)G (] 

•/ —OO 



(Re 



+21, • I, fdt [Gq(R, e) + Gi(R, e)] 

«/ —OO 



(39) 



Here, the trace over spin degrees of freedom were calcu- 
lated by using the following identities 

Tr{(I,. S )(I,. S )} = 2I l .I„ (40) 
Tr{(I, • 5)5,(1, • «)} = -2i(Ii x I,-),, (41) 
Tr{(I, • S ) Sz (I, • S )s z } = 2(2/? 7/ - I,: • I,). (42) 

All integrals in Eq. (|39| are of the same type, 



Im 



de 



p ikr,\z\ ik t\z\ hiTrr 

, V n+ * . n+ = ^siikrnM), (43) 
k n + l0 + Krj' + iv + kc 
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and can be easily evaluated by changing variables from 
the original e to k„ + k„i . We denote the real part of 
the sum of two Fermi wavevectors as k VtV > = Re[k v (eF) + 
fc^'(ef)] with the indices 77,77' = ±1. As a result, we 
arrive at the RKKY Hamiltonian in the form of Eq. (|2|) , 
where the components of the spin susceptibility tensor 

ij 
a/3 



= Xa@ (R) are explicitly given by 



2irhvF 



Xxx — 



2nhvF 



si(fc+ 1+ M)+si(fc-,-|z|) 
4-2cos(2K-R)si(fc_ i+ |z|) 
2si(fc_. + |z|) + cos(2K-R) 



(44) 



x [si(fc+, + |z|) + si(/c_,_|z|)]), (45) 



Xxy 



a 2 k G 



2nhvF 



sin(2K • R) 

x si(A + , + |z|)-si(fc_,_|z|)l, (46) 



with Xxx = Xyy, Xxy = ~Xyx, and all other compo- 
nents being zero. First, we note that the off-diagonal 
components Xxy,Xyx are non-zero. They describe the 
response to a perturbation applied perpendicular to the 
nanotube axis. This opens up the possibility to test the 
presence of SOI in the system by measuring off-diagonal 
components of the spin susceptibility tensor Xa/3- Sec- 
ond, the spin response in a direction perpendicular to 
the z-axis cannot be caused by a perturbation along the 
z-axis, thus Xyz = Xxz = 0. This simply reflects the 
rotation-invariance of CNTs around their axes. The dif- 
ference between the diagonal elements of the spin suscep- 
tibility tensor, Xxx = Xyy and Xzz, again arises from the 
SOI and is another manifestation of the broken rotation 
invariance of spin space. In total this means that the 
RKKY interaction given in Eq. ([39]) is anisotropic in the 
presence of SOI, giving rise to an Ising term oc and 
a Moryia-Dzyaloshinsky term oc (Ij x lj) z , in addition to 
the isotropic Heisenbcrg term oc Ij ■ Ij . 

Quite remarkably, when the Fermi level is tuned close 
to the gap opened by the SOI, then the off-diagonal and 
diagonal components of the susceptibility tensor become 
of comparable magnitude, see Figs. [5] and [3] This has im- 
portant consequences for a Kondo lattice system, where a 
highly anisotropic RKKY interaction will give rise to an 
ordered magnetic phase even at finite temperatures i21 As 
a potential candidate for such a Kondo lattice^ we might 
mention a CNT made out of the 13 C-isotope^ where 
each site of the graphene lattice contains a nuclear spin- 
1/2 to which the itinerant electrons couple via hypcrfine 
interaction^ 

We note that the susceptibility depends on two Fermi 
wavevectors via k v ,y in a rather complicated way (see 
Figs. Hand El). In the absence of SOI, we recover the 
result for the spin susceptibility on the same sublattice, 
X AA (£V) [see Eq. ([2"3])]. The leading term in the spin 
susceptibility for different sublatticcs, x j4S (R) 5 can be 




Xxy/XO 



500 

z (nm) 



FIG. 3. The same as Fig.[2](a,b) but with the chemical poten- 
tial (j, = 474 meV being tuned above the gap opened by SOI. 
The spin susceptibility decays as 1/z and exhibits beatings, 
with period determined by the SOI parameters. 



obtained from Eqs. (|44[) - (|46[) by putting a minus sign 
in front of cos(2K • R). In addition, as shown in Sec. 
IIV B[ the spin susceptibility vanishes if the Fermi level 
is tuned inside the gap in semiconducting CNTs, so that 
both k±{ep) are purely imaginary. If the chemical po- 
tential is inside the gap opened by the SOI, the Fermi 
wavevector fc_(e^) is still imaginary, at the same time 
k + {tp) is real, giving fc +:+ = 2k + {ep), = 0, and 

= k + (ep). This results in the behavior of the 
spin susceptibility shown in Fig. [3J The strength of the 
RKKY interaction decays oscillating as 1/R. The oscilla- 
tion period is determined by fe+ + for Xxy and by fc+,- for 
Xxx and Xzz , see Fig. H If the chemical potential is above 
the SOI gap, then both wavevectors k±(ep) are real, giv- 
ing rise to oscillations with two different frequencies that 
result in beating patterns for the spin susceptibility, see 
Fig. El 



VI. RKKY WITH MAGNETIC FIELD 

In Sec. |V]wc demonstrated that the presence of SOI, 
which breaks rotation invariance in spin space, leads to 
an anisotropic spin susceptibility. Another way to break 
this rotation invariance is to apply a magnetic field, which 
also breaks time-reversal invariance. In this section we 
again neglect sublattice asymmetries discussed above and 
focus on the effects of a magnetic field B applied along 
the nanotube axis z for a semiconducting nanotube. The 
Zeeman term Hz = AzSz = HbBs z lifts the spin de- 
generacy. Here, /is is the Bohr magneton. The orbital 
term H or b = "fhvpk ma g leads to a shift of the trans- 
verse wavevector k G by k mag , which finds its origin in the 
Aharonov-Bohm effect and is given by k mag = TrBd\e\/2h 
for a nanotube of diameter d. Thus, the valley degen- 
eracy of the levels is also lifted. The spectrum of the 
effective Hamiltonian Ho + Hz + H or i, in the case of a 
semiconducting nanotube is given by 

e 7iS = hv F k Gn + s( 7 /3+ + A 2 ) + Hv F k 2 /2k G . /y , (47) 
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where fc(j. 7 = k G + jk mag . The Green functions in mo- 
mentum space can be found similar to Eq. (|32[) . As a 
result, we arrive at the following expression for the Green 
functions, 

G(k, e + i0 + ) = ^[G (7, k, e) + Gi( 7 , k, e)s z ], (48) 

7 

where 

kg,-/ 



Gi(7,k,e) 



E 



flVF 

revec 

energy e from Eq. (|47|) as 



E 



"7,s 



(49) 
(50) 



Here, we define wavevectors function of the 



' 2k Gn [e - hv F k Gn - s( 7 /3+ + A,)] 



Hvf 



(51) 



which can take both non-negative real and imaginary 
values. The Green functions in real space are found by 
Fourier transformation, 



G (R,e) 
Gi(R,e) 



V- afc G,7 c »7K-R e "^ N 

^ 2hv F 

7.s 



7,s 



i0+ 



E 

7,s 



7 e i 7 K-R 



SC. 



z0+ 



(52) 
(53) 



By substituting Eqs. ([52]) and ((53]) into Eq. ([39]) . 
we arrive at the effective RKKY Hamiltonian. Since 
kmag/kc « 1, we can neglect the dependence of the 
spectrum slope on the magnetic field, which simplifies 
the calculations considerably. 

At the end we arrive at the following expressions for 
the spin susceptibility tensor components 



Xxx 



H4 



E [ si ( fc -> 



\z\) + cos(2K • R)si(fc 7 



7,s 



X 



a 2 k G 



A-kHvf 

X ^ ^ ^Si(fc-y jS ;-y jS | 



:|) + cos[2K • R]si(fc 7jS;7 



7,s 



Xxy 



a 2 k G 



V7ssin[2K • R]si(fc 7 



z\) , (54) 



(55) 
(56) 



Xyy, Xxy = ~Xyx, and the rest being zero. 



with Xx 

Here, we use the notation /c 7 . s;7 '. S ' = Re[>t 7 . s (ei?) + 
k 7 ' Again, the RKKY interaction decays at large 
distances as 1/R. The spin susceptibility also exhibits 
oscillations and beating patterns (similar to ones shown 
in Fig. |3]) determined by four different Fermi wavevec- 
tors k ls (eF), see Fig. 2J Finally, we note that the same 



0.4 
0.2 



.... Xxy /Xo 




FIG. 4. The same as Fig. [2] (a,b) but with the chemical po- 
tential /j, = 472.7 meV being tuned in such a way that there 
are three pairs of states at the Fermi level and with a mag- 
netic field B = 1 T. The spin susceptibility decays as 1/z 
and exhibits beatings, with period determined by the SOI 
parameters and the magnetic field. 



beating patterns arises also for nanowires with parabolic 
spectrum in the presence of a magnetic field giving rise 
to a Zeeman splitting. 



VII. SPIN FLUCTUATIONS 



The spin susceptibility is a fundamental character- 
istics of the system. At zero frequency it describes 
the RKKY interaction between localized spins. At fi- 
nite frequencies oj the spin susceptibility gives access to 
the equilibrium spin noise in the system. For a gen- 
eral observable, s a , the fluctuation-dissipation theorem^ 
connects the dynamical structure factor S aa (q,uj) = 
(l/2ir)fdt (so,(q, t)s a (— q, 0))e Mt , describing equilib- 
rium fluctuations, with the linear response susceptibility, 
Xaa (q.,oj) = (-i/H)fdt e(t)([« a (q,t),« a (-q,0)])e &rt , 



Imx QQ (q,w) = -(e k s T - l)5 QQ (q,a;) 



(57) 



for u> > 0. Below we calculate the spin susceptibility 
for both metallic and semiconducting naiiotubes at finite 
frequencies, and obtain via Eq. ([57]) the spin correla- 
tion function (spin noise) S' aQ (R, uj) in real space at zero 
temperature, 



S QQ (R,w) = Imx Qa (R,w) 

7T 



(58) 



Here, we assumed that the system is invariant under par- 
ity transformation so that x(Qi w ) becomes an even func- 
tion of q. For simplicity, we consider only the case with- 
out SOL 
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A. Metallic nanotubes 

The spin susceptibility at finite Matsubara frequencies 
uj„ is given by 

2 r°° 

Xo(R, iwn) = l J dr G (R, r)G (-R, -r)e^ T , 

(59) 

where we modified Eq. ((5|) accordingly, and the Green 
function Gq(R, t) was found before [see Eqs. ([To]) and 
(|T7|) ]. For the spin susceptibility on the same sublatticc 
we get 

AA/t> ... . \ 2a 2 2/ 



X^(R,ia; ri ) = ^-cos 2 (K.R) 



: / dr e 1 ' 
lo 



T [z sin(fcFz)] 2 — [vft cos(kFz)]' 



Introducing the notation 



I n {x) 



[(vft) 2 + z 



dr 



212 



(60) 



lo (1 + r 2 ) 2 ' 
where n = 0, 2, the susceptibility can be rewritten as 
2a 2 



(61) 



X^(R,iW n ) 



sin 2 (fcf z)/q 



TT 2 Hvf\z 

w n |z| 



cos 2 (K-R) 



cos 2 (kpz)l2 



u n \z\ 

Vf 



(62) 



The asymptotics of I n (x) for small x is given by Iq(x) « 
7r/4 + £x/2 and Iz(x) ps tt(1 — 2x)/4 — £x(lnx + 7 e — 
1/2), where 7 e is the Euler constant. Performing the 
standard analytic continuation from Matsubara to real 
frequencies^ we obtain the spin susceptibility Xo A { U) )i 
and from its imaginary part the spin noise, see Eq. (|58|) . 
Explicitly, the dynamical structure factor at zero tem- 
perature and for vf 3> ui\z\ > is given by 



2a 2 u 



„ cos 2 (K • R) cos 2 (fc F z). (63) 



For the susceptibility on different sublatticcs we find 
2a 2 



TT 2 hv 



■sin z (K • R) 



\ T / w„|z|\ 2 / w„|z| 

sin (kFZjh — — 1 — cos (k F z)I ' 



vf 



Vf 



(64) 



with the corresponding dynamical structure factor being 
given by 



S AB (K,cu) 



2a 2 w 
ir 2 Vp 



sin 2 (K-R) sm 2 (k F z). (65) 



The dynamical structure factors S AA (R,u)) and 
S AB (R 7 uj) are linear in u>, and, as expected, vanish at 



zero frequency (we recall that we work at zero temper- 
ature). Moreover, the noise is strongly suppressed at 
some special points R z in space that satisfy the condi- 
tion kpz = mr/2, where n is an integer. 

In the opposite limit ui\z\ 3> vp we use the follow- 
ing asymptotics, Iq(x) « i/x + irxe~ x /2 and ^2(2;) ~ 
—2i/x 3 — irxe~ x /2 for x ^ 1. After analytic continua- 
tion, the dynamical structure factor is then given by 



S AA (R, W ) 



S AB (R,w) 



ir 2 v 2 F 



\ 2 v 2 F 



cos 2 (K-R)cos( ^ 

Vp 

sin"(K-R)cos( ^ 
1 vf 



(66) 
(67) 



In this regime, S AA (R,,ui) and S'" 4B (R, w) arc not only 
linearly proportional to the frequency but also oscillate 
rapidly as a function of frequency. This implies that in 
real time the spin noise is only non-zero for times t and 
distances z satisfying \z\ = Dpt. 



B. Semiconducting nanotubes 

For semiconducting CNTs all calculations for the fre- 
quency dependent susceptibility arc similar to the ones 
for onc-dimcnsional systems with parabolic spectrum, 
being available in the literature^ At zero temperature 
ImXaa(<7,w) is given by 

k 

X 5(hw +£fc,s,7 ~ 6fe+5,*,7')- (68) 

The upper and lower frequencies lj± arc defined as 
VFq 2 vpkF 



w±(g) 



± 



2k G k G 



(69) 



The imaginary part of the spin susceptibility is non-zero 
only for frequencies 



< M < uJ+{q) 

and is given by 

ImXaa(g,7)7') w ) 



ak G 



4lHvf\<1 



(70) 



sgn(w). (71) 



To arrive at the expression in real space we perform the 
Fourier transformation, 

Imx aa (R, oj) = (1 ± cos(2K • R)) 

where the range of the ^-integration is determined from 
Eq. (|70[) . and the positive (negative) sign corresponds to 
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Xaa (Xaa)- For high frequencies, w > v F k 2 F /2k G , the 
dynamical structure factor is given by 

£(R,w) = (l±cos(2K.R)) ; 



2tt 2 vf 

x [sKq+Hlzl) - si(«_(w)|«|)] sgn(^), (73) 

where wavevectors q±(ui) are positive solutions of the 
equations |cj| = |w±(g±)|. For low frequencies, < uj < 



VFk F /2kc, the same equation 



;_ (3- i)| has three 



non- negative solutions,— q—,\ < Q— ,2 < q~. 3- In this case 
5(R, w) is given by 



5(R,w) = (l±cos(2K-R))- 



-sgn(w) 



2tt 2 vf 

x [si(g + (w)|z|)-si(g_,i(o;)|z|) 

+ a(q-,2(u)\z\)-a(q-,3(u)\z\)]. (74) 

We note that the expression is composed of several con- 
tributions and thus leads to beating patterns of the spin 
noise, similar to the one before for the spin susceptibility. 



VIII. GRAPHENE NANORIBBONS 

A. The effective Hamiltonian 

In the last part of this work, we turn to graphenc 
nanoribbons, which arc finite-size sheets of graphenc: 44 
The nanoribbon is assumed to be aligned along the z- 
direction and to have a finite width W = Na in IE- 
direction, with N being the number of unit cells in this 
transverse direction. Here, we focus on armchair nanorib- 
bons, characterized by the fact that the x-axis points 
along one of the translation vectors of the graphene lat- 
tice. The effective Hamiltonian is given by 



H = hv F {jk x ai + k z a 2 ) 



(75) 



which determines the low-energy spectrum around the 
two Dirac points K = — K' = (47r/3a, 0). Here, k z is 
the momentum in z-direction. The momentum k x in 
x-direction is quantized due to the vanishing boundary 
conditions imposed on the extended nanoribbon4^ If the 
width of the GNR is such that N = 3AI + 1, where M 
is a positive integer, the GNR is metallic with k x = 0. 
Otherwise, the nanoribbon is semiconducting. 

The eigenstates arc written as ip = ~}2 ai 4>a-ye l " ,K ^ x : 

* = {4>ak, 4>bkAak> ,<t>BK'), where a = A,B. The cor- 
responding spectrum and wavefunctions that satisfy the 
vanishing boundary conditions (for ip) are given by 



e,k z 



1 k - ; 



HCMC-i), 



e c = (hv F k z 
for a metallic GNR, and 



(76) 
(77) 



e.k z 



e ik z z ^ e i<p s +ixk™ n ^ix 
^itps—ixk^ 



e± = ±hv F y / [k 



rain \ 2 



kl 



cfcr "), (78) 
(79) 



for a semiconducting GNR. Here, C, = ±1 is the eigen- 
value of the Pauli matrix 02, and we use the nota- 
tion e iv ° = [k™" 1 - ik z ]/y/{k™ in ) 2 + k 2 z , with = 
7r/3(iV + 2)a. 



B. Spin susceptibility 

1. Without SOI 

To calculate the spin susceptibility for a metallic 
nanoribbon that has a linear spectrum given by Eq. (|77l) , 
we again work in the imaginary time formalism, see Eq. 
([5])- The calculations are quite similar to the ones pre- 
sented before in Sec. IIV Al The only change in the ex- 
pressions for the spin susceptibility in comparison with a 
metallic nanotubc [see Eqs. (fT5|) and (|19|) ] is in the fast 
oscillating prcfactor, 



AB, 



2TThvp\z\ 



x sin 2 (K • Ri) sin 2 (K • R,-) cos(2fc F z). (80) 

Similarly, for semiconducting nanoribbons the spin sus- 
ceptibility is given by Eqs. (|24|) and (f27|) . where the 
fast oscillating prcfactors 1 ± cos(2K • R) are replaced by 
sin^K-ROsin^K-Rj). 



2. With SOI 

The intrinsic SOI in graphenc is only several fieV, so 
it is rather weak. Moreover, the Rashba SOI generated 
by an externally applied electric field E is in the range 
of tenths of [iey for E = 1 V/nm.^ Such small SOI 
values might be hard to observe. However, the Rashba 
SOI generated by a spatially varying magnetic field opens 
new perspectives for spintronics in graphene.— In this 
case, the SOI strength can be exceptionally large, reach- 
ing hundreds of meV. A nanoribbon in the presence of a 
rotating magnetic field with period 27r/fc„, described by 
the Zccman Hamiltonian 



Hi 



A z [s y cos(k n z) + s z sin(fc„z)] , (81) 



is equivalent to a nanoribbon with Rashba SOI in the 
presence of a uniform magnetic field, 

H 1 ^ = Ul(H + H^)U n = H + A z s y + A n so s x a 2 , (82) 

where the unitary gauge transformation is given by U n = 
cxp(ik n zs x /2). The period of the magnetic field deter- 
mines the strength of the Rashba SOI A™ = hvpk n /2, 
while the amplitudes of the uniform and the rotating 
fields are the same and given by A^. 10 ' 13 The spectrum 
of a metallic (semiconducting) GNR in the presence of 
such SOI consists of two cones (parabolas) shifted along 
the momentum axis against each other by k n /2. Ev- 
ery branch of the spectrum possesses a well-defined spin 
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polarization perpendicular to the z-axis that is along the 
nanoribbon. A uniform magnetic field only slightly modi- 
fies the spectrum by opening a gap at zero momentum. In 
the following discussion we neglect the uniform magnetic 
field working in the regime where the induced Rashba 
SOI is stronger than the Zeeman energy 2A^. 

As a result, similar to the semiconducting nanowirc, 
one can gauge away the momentum shifts by rotating 
the spin coordinate system as follows, 

I X (R) = I x cos(fc„z) + I y sin(fc„z), (83) 
I y (B.) = I y cos(fc„z) - I x sin(fc„z), (84) 
I»(R)=i». (85) 

The same transformation should be applied to the elec- 
tron spin operators hs/2. The effective RKKY Hamilto- 
nian in this rotated coordinate system is the same as in 
the system without SOI and is given by Eq. (|80l) . To 
return to the laboratory frame, we perform the following 
change 

Ii • Ij -> Ii(Rj) • Ij(Rj) = cos(fc„z)Ii • Ij 

+ [1 - cos{k n z)]I?I] - sm{k n z)(Ii x lj) z . (86) 

The spin susceptibility tensor has non-vanishing off- 
diagonal components, which, again, indicate a broken in- 
variancc of spin space induced by the magnetic field or 
the Rashba SOI. As before, this gives rise to anisotropic 
RKKY interactions of Ising and Moryia-Dzyaloshinski 
form. 



IX. CONCLUSIONS 

In the present work we studied the Rudcrmann-Kittcl- 
Kasuya-Yosida (RKKY) interaction in carbon nanotubes 
and graphene nanoribbons at zero temperature in the 
presence of spin orbit interaction. Our main results are 
summarized in the following. 

The spin susceptibility in metallic CNTs, characterized 
by a Dirac spectrum (gapless and linear), crucially de- 
pends on whether the localized spins that interact with 
each other are from the same or from different sublat- 
tices. In particular, if the Fermi level is tuned exactly to 
the Dirac point where the chemical potential is zero the 
interaction is of ferromagnetic type for spins on A- A or 
B-B lattice sites, whereas it is of antifcrromagnctic type 
for spins on A-B lattice sites. In semiconducting CNTs, 
with a sizable bandgap, the spin susceptibility depends 
only slightly on the sublattices. In all cases, the spin sus- 
ceptibility is an oscillating function that decays as l/R, 
where R is the distance between the localized spins. 



The spin orbit interaction breaks the spin degeneracy 
of the spectrum and the direction invariance of the spin 
space. As a result, the spin susceptibility is described 
by the tensor Xap that has two non-zero off-diagonal 
components \xy = ~Xyx, the finite values of which sig- 
nal the presence of SOI in the system. Moreover, the 
RKKY interaction is also anisotropic in the diagonal 
terms, x*z ^ Xxx = Xyy Quite surprisingly, we find 
that all non-zero components, diagonal and off-diagonal, 
can be tuned to be of equal strength by adjusting the 
Fermi level. These anisotropies, giving rise to Ising and 
Moriya-Dzyaloshinski RKKY interactions, thus open the 
possibility to have magnetic order in low-dimensional sys- 
tems at finite temperature^ 

We note that, in contrast to semiconducting nanowires, 
the SOI cannot be gauged away by a unitary transforma- 
tion in CNTs, giving rise to a more complicated depen- 
dence of Xufi on the SOI parameters. In the same way, a 
magnetic field along the CNT axis breaks both the spin 
and the valley degeneracy, leading to a dependence of the 
spin susceptibility on four different Fermi wavevectors. 

The spin susceptibility at finite frequencies also al- 
lows us to analyze the spin noise in the system via the 
fluctuation-dissipation theorem. We find that the dy- 
namical structure factor S aa (R, uj) is linear in frequency 
and oscillates in real space. 

Metallic armchair GNRs behave similarly to metallic 
CNTs. Indeed, in both cases the spin susceptibility shows 
a strong dependence on the sublattices, with, however, 
different fast oscillating prcfactors. A Rashba-like SOI 
interaction can be generated in armchair GNR by peri- 
odic magnetic fields. In contrast to CNTs with intrinsic 
SOI, this field-generated SOI can be gauged away giv- 
ing rise to a simple structure of the spin susceptibility 
tensor4£ 

In this work we have ignored interaction effects. How- 
ever, it is well-known that in one- and two-dimensional 
systems electron-electron interactions can lead to inter- 
esting modifications of the spin susceptibility, for in- 
stance with a slower power law decay such as l/R 9 , with 
< g < 1 in a Luttingcr liquid approach to interact- 
ing one-dimensional wires£ii£ It would be interesting to 
extend the present analysis and to allow for interaction 
effects^ in the spin susceptibility for carbon based mate- 
rials in the presence of spin orbit interaction, in particular 
for metallic CNTs and GNRs at the Dirac point. 
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